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Abstract. Exterior calculus and moving frames are used to describe curved elastic 
shells. The kinematics follow from the Lie-derivative on forms whereas the dynamics 
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C^** ' via stress-forms. 
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O ■ 1- Introduction 

^ ' Shell theory reduces the three-dimensional elasticity of a shell to an effective two- 
■*—> 

dimensional theory of its middle-section. This article presents a new short derivation of 
the shell equations of Kirchoff-Love based on the method of moving frames and exterior 
^ ' calculus. 



Early shell theory dealt with church bells. Lord Rayleigh studied vibrating shells 
motivated by the question whether or not a full musical scale could be played on bells pp. 
Vibrations of shells for airplanes and cars continue to be important in civil engineering. 
There, the numerical modelling of shells remain an active area of research relevant for 
e.g. car crashes. Also computer scientists have joined the efforts in order to present 
realistic computer graphics [2]. Lately, shell theory has been applied to cell membranes 
[Sj and even nanotubes j3j. 

There are several approaches to shell theory. Kraus presents a variational method 
in Equivalently, assuming Newton's second law for an infinitesimal volume element 
and integrating over the thickness of the shell lead to the same results, see e.g. jE]. 
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Unfortunately the calculations are complicated: the dynamic part of this well-known law 
involves a covariant derivative of the stress tensor and requires intricate manipulations 
of Christoffel symbols for the full shell, respectively its middle-section. 

However, a general observation is that complexities in tensor calculus can be avoided 
if it is possible to work in an orthonormal frame, i.e. in the case of Riemannian 
manifolds. Then the number of connection coefficients are drastically reduced and 
calculations become simpler and more transparent. For example, Maxwell's equations in 
electrodynamics and calculations of curvature in general relativity become much simpler 
when using the exterior calculus of differential forms and Cartan's method of moving 
frames [Zj. It is in this spirit we shall discuss shell theory. 

1.1. Thin shell assumptions 

For simplicity, we follow the assumptions of Kirchoff and Love for a thin shell: vertical 
coordinate z and intrinsic coordinates a±, a 2 : 

(i) thickness small compared to radii of curvature 

(ii) small displacements 

(iii) vanishing normal stress 033 = 

(iv) preservation of normals e i3 = 

(v) linear dependence of membrane field U l = u l (a l , a 2 , t) + z/3 l (a l , a 2 , t) 

(vi) constant dependence of flexural field W = ^(a 1 , a 2 , t) 

1.2. Outline 

The presentation centers around a proof of the shell equations. Some familiarity with 
differential forms and continuum mechanics is assumed. Sect. El introduces the geometry 
via differential forms, sect. El treats the kinematics of shell deformation and sect. 0|the 
equations of motion. A summary and some concluding remarks will be given in the last 
section. 
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2. Geometry 

2.1. Geometry of the mid-section AA 

The shell equations must be covariant and therefore it is enough to validate them in 
one set of coordinates. Hence, on the mid-section AA we can choose lines of curvature 
coordinates a 1 and a 2 having axes aligned with the principal directions of curvature. 
Such coordinates can always be chosen locally on a two-dimensional surface |Bj. When 
the radii of curvature differ Ri ^ R 2 the axes can be chosen uniquely. Otherwise, at 
an umbilic point R\ = R 2 ^ oo or at a planar point R\ = R 2 = oo, we shall assume a 
choice made. Write in these coordinates, the mid-section metric in the Lame-form 

ds 2 = A\ (da 1 ) 2 + A\ (da 2 ) 2 . (1) 

The Ai are called Lame-parameters. Instead of working directly with the coordinates a 1 




Figure 1. Middle section segment with curvature radii i?i,i?2 and increments 
A 1 da 1 ,A 2 da 2 . 



and a 2 as done traditionally we shall work with a frame derived from these particular 
coordinates. Thus define an orthonormal frame on the mid-section via the forms 

4> a = A a da a - (2) 
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Then the metric of the mid-section the first fundamental tensor takes the simple form 

a = g ab dx a <g> dx b = cf) 1 ® <p l + 4> 2 ® <p 2 = 5 ab (p a ® b (3) 

and the curvature tensor of the mid-section the second fundamental tensor becomes 

d = d ab dx a ®dx b = ^0 1 ®0 1 + ^-0 2 ®^ 2 . (4) 
H\ R 2 

£.2. Geometry of the shell S 

Denote the coordinates on the three-dimensional shell S: a 1 , a 2 and z. The latter 
coordinate is perpendicular to the middle section of the shell, see fig. [0 The direction 
of the z-axis follows that of j3] but is opposite from e.g. [UJ E]. The thickness of the 
shell is called h. z equal ±/i/2 defines the upper respective the lower section of the shell. 
Likewise z = corresponds to the mid-section A4. The shell metric is 

ds 2 = A\(l + z/R^ida 1 ) 2 + A\{\ + z/R 2 ) 2 (da 2 ) 2 + dz 2 . (5) 

The factors 1 + z/Ri are understood as follows: An increment on the middle surface in 
the direction i spans the angle 

dtp 1 = Aidot jRi 
so the increment at the elevated position is 

dx\z) = (Ri + z) dtp* = Ai(l + z/R i )da i . 
Finally, we shall assume furthermore the existence of derivatives of Ai, Rj. 

2.2.1. Orthornormal frame of S We introduce a three-dimensional frame as done for 
the mid-section instead of working with the coordinates a 1 , a 2 and z as in . Thus, 

from © define an orthonormal basis of one-forms as: 

1 = A x (1 + z/R x ) da 1 , 2 = A 2 (1 + z/R 2 ) da 2 and 3 = dz . (6) 
From the one-forms construct the volume element 

dV = 1 A 2 A 3 = A X A 2 (1 + zjR x ){\ + z/R 2 ) da 1 A da 2 A dz (7) 
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The Hodge-dual forms to the one- forms are oriented area-elements 

dS a = *0 a (8) 

and fulfill 

@ a A dS 6 = 5 a b dV . (9) 

In component form 

dS, = e ijk @ j ® & k (10) 

with Eijk the totally antisymmetric tensor. Note, in the following we reserve for 
another tensor, the strain, and hopefully no confusion will arise. As the area-forms are 
not closed cZ(dSj) ^ they are not exact dS.; ^ d(l — form), and hence the notation of 
these area forms is slightly misleading although conventional. 

2.2.2. The connection of S The one-forms of the frame lead to the connection 
coefficients via the condition of vanishing torsion [7j: 

d& + <J j A @ j = . (11) 

Calculating d® 1 from © the connection form ui l j is found by inspection: 

u = ( W y (12) 

n (A 1 (l+z/R 1 )), 2 & 1 -(A 2 (l+z/R 2 )) il & 2 1 v 

U A 1 A 2 {l+z/R 1 )(l+z/R 2 ) Ri+z \ 

-(A 1 (i+^/fli)),2e 1 +(A 2 (i+ z /fl 2 )) ;1 e 2 g e 2 



A 1 A 2 (l+z/R 1 )(l+z/R 2 ) " R 2 +z 

V — ®i- ®i- J 

Ri+z R 2 +z 

/ n A 1 , 2 ® 1 /{l+z/R 1 )-A 2 , 1 ® 2 /{l+z/R 2 ) , 

/ U A\A 2 R x +z \ 

-A 1 , 2 Q 1 /(i+ z /fli)+A 2 ,ie 2 /(i+ z /R2) q e 2 



A±A 2 R 2 +z 

V — ®i- o / 

where derivatives are denoted with a comma Of /da 1 = fj. The final simplification in 

()12|) follows from the Gauss-Codazzi equations, shortly discussed in the following. 

Finally, note that in the subsequent calculation of various covariant derivatives the 
number of connection coefficients T l k - = L ek uj l j are far less in an orthonormal frame as 
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compared to a general coordinate basis as these coefficients for one pair of indicices now 
are antisymmetric and not just symmetric. 



2.2.3. Gauss- Codazzi There exist constraints on the functions At and Ri in the metric 
(jSJ for defining a valid surface. These equations come about by expressing the curvature 
via the connection. Even though the shell is curved the metric (JSJ) is just a metric for 
the flat three dimensional space 1R 3 . Consequently the curvature is zero: 

= nv = dJj + w\ A u k j (13) 

For instance for % = 1 and j = 3 

du} x z = (A 1 /R 1 ) i2 da 2 A da 1 (14) 



and 



u \ A u\ = (Al(1 + " /jRl)) ' 2 da' A da 2 , (15) 
i?2 + Z 



so 



;; ' 2 -'^i (l + z/R 2 ). (16) 



At the middle surface z = the classical Codazzi equation is found 
A h2 Mi 



i?2 \Rl 
Applying (|17jl to (fTHj) gives 



;i7) 



+ z/i?i)), 2 = A h2 (1 + ^/,R 2 ) (18) 

used in (II 2j) . Likewise other identities between Ai, A 2 ,Ri and _R 2 hold generalizing the 
Gauss-Codazzi-equations to the shell metric. By considering i — 1 and j = 2 the result 
of Gauss follows: 

Hence using (ITHjl . the total curvature = related to the embedding may be 

expressed entirely in terms of the intrinsic functions Ai belonging to metric of the 
surface. This is the lines of curvature coordinates version of Gauss's theorema egregium 
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2.2.4- Covariant Derivative Having denned the connection u>, the covariant derivative 
acting on forms is 

V = d + uj, (20) 

where d is the flat exterior derivative and u> acts from the left with the wedge-product. 
From the condition of no torsion (JTTj) follows that the covariant derivative vanish 
on frame co-vectors respective tangent-vectors. The latter holds since the metric is 
covariantly constant. Thus 

ve* = o. (21) 

To exploit this structure we shall consequently define quantities of interest as forms. 

2.2.5. The derived connection of Ai The mid-section Ai is embedded in the shell S: 
i:M^S (22) 

so restricting differential forms to the mid-section is the pull-back 

i* : A*(S) — > A*(M) . (23) 

This restriction is obtained by putting z — 0. For instance 

i*(0 a ) = a and i*0 3 = O. (24) 

Likewise, pulling back the connection of the three-dimensional shell u) l j to the mid- 
section i* gives the two-dimensional connection u>% of Ai. Namely 

d(b a = -u\\ z=0 A(f> b = -Z a b A cp b (25) 

with a,b = 1,2 only, where the "tilde" shall be used in the following when referring to 
the mid-section. Formally f)25|) holds since the exterior derivative commutes with the 
pull-back 

di* = i*d (26) 

and 

d<f> a = d{t*& a ) = t*d@ a = fi-oJ^ A &) = -t*{u: a b ) A (j) b . (27) 
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Thus, all connection coefficients of the mid-section T l - k are contained in the single one- 



form 




A^tft 1 - A 2il <f> 2 
A±A 2 



(28) 



where 



~a pa ic — b 



(29) 



Finally, these connection coefficients allow us to define the covariant derivative on the 
mid-section, V. 

3. Shell kinematics 

This section describes the kinematics of a shell. The appropriate deformation field is 
the strain field e. This field corresponds to relative displacements or extensions. 

3.1. The displacement field U 

The displacement is a tangent vector field and hence contra- variant. According to 
Kirchoff-Love: 



More general functional forms exist: e.g. quadratic dependence in the flexural field w 
and even general Taylor expansions in z have been proposed by Vlasov j^J. 

3.1.1. Rotation angles f3 l If the (3 in (jHUj) is sufficiently small its components play the 
role of angles. Thus when seen from the middle section the displacement at height z in 
the direction % equals zdip 1 ~ z(3 l . However, in Kirchoff-Love's theory, the rotations f3 l 
are not considered as independent fields but can be found from u % and the gradient of 
w. This we show in the following sect. 13.31 using the assumption of vanishing normal 
strain = 0. 



U = U 1 e 1 + U 2 e 2 + We 3 = (u 1 + z^)ex + (u 2 + zf3 2 )e 2 + we 3 



(30) 
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3.2. Strain field: connection to the Lie- derivative 

In linear elasticity the strain field follows from the displacement field by an application 
of a Lie- derivative on the metric 

e = e ij @ i ®@ j = ^£u(g) with g = gijdx* <g> dx j = 5 jk @ j <8> @ k (31) 

The Lie-derivative £ is with respect to the physical displacement field U. 

The strain field is calculated with respect to the orthonormal frame (JBJ). As £ is a 
derivation: 

£(& ® @ k ) = £{@ j ) ® @ k + e J ® £(0 fe ) (32) 

it suffices to investigate its action on a co- vector @ J . Using Cartan's formula JT] for 
the Lie-derivative and ([2 lft: 

£u(0 J ) = (d iu + tu d)0 j = dC/ J ' + Lu{d@ j ) = dU j - L V (u) j k A @ k ) 

= dU j + U k ^ k - {iv^ k )Q> k = VW - (kju^)©" • (33) 
The antisymmetry uj^ k = —u) k j gives the familiar result from continuum mechanics: 

e = -(W ® e fc + e fc ® vu j )g jk = -(v <g> e^- + t/ j e- ? ' ® v) = -(vu + uv) .(34) 

In (}3"4*j) the metric lowers the indices of the displacement field coordinates from a tangent 
vector to a co- vector. Sometimes this is made explicit by writing U = U^. 

3.3. Lie- derivative gives variation of curvature 

We wish to emphasize the connection of the strain in shells to the variation of the second 
fundamental tensor. This is done using the Lie-derivative instead of the conventional 
expression of the strain (}3"4*j) . 

Expanding the shell metric To find the strain with respect to the mid-section the full 
shell metric (JSJ) is written out in terms of simpler mid-section one-forms 

g = (f) 1 ® (f) 1 + 2 ® d> 2 + 3 ® 3 

+ 2z (Ltf^tf + Ltf®^ +z z (L^^^ + ^tf®^ 

= a + e 3 ® e 3 + 2^d + z 2 d 2 (35) 
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with d 2 b = d ac dl. Deriving with the Lie- derivative gives 

e =2 £ u(g) 

= - £u(a) + - (Vw <g> 3 + 3 ® Vic) + Cv(z d) + - £u(^ 2 d 2 ) 
= - £u(a) +wd + z £u(d) + z w d 2 + i (Vto ® 3 + 3 ® Vto) + 0(z 2 ) 

(36) 

Above, the Lie-derivative contracts the mid-section forms a via the interior derivative 
ljj. Hence, it is advantageous to express the frame coordinates of the displacement field 
with respect to the mid-section: 



[l + z/Ri)Ai da 1 
= (u' + z mi -zdl)^^- + 0(z 2 ) 

= (u* + z{fi - d) u 1 )) e~i + 0{z 2 ) = (u l + zp l )e t + 0(z 2 ) (37) 

with 

u = u (38) 

/3 = /3-d-u. (39) 
To proceed consider the first term in ()36|) 

i £u(a) = l - £ u+z/3 (a) = i £ u (a) + i £ 2/3 (a) (40) 
The first term in (J4*U|) reduces to a membrane strain expressed in mid-section coordinates 

^ £u(a) = ^(u a]b + u b]a ) (f> a <g> </> b = ^(V <g> u + u ® V) (41) 

with the corresponding covariant derivative. The second term in (|4U|) gives a normal 
strain plus a strain similar to ()41|) of order O(z): 

^£ 2/3 (a) = ^ a (^®0 b + 6 ®^)cU + ^(V® / 3 + / 3® V) 

= i(0 3 ®/3 + /3®0 3 ) + z^(V®/3 + /3® V) (42) 

_ _ 

As |jj = the normal strain €33 vanishes in (|3fi|) . So effectively 

Vw = = dw(a , a 2 ) = Vw . (43) 
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Then the requirement of vanishing normal shear strain e a3 = with (|36|) and (|42j) 
reduces to 

/3 = -Vu> + 0{z) . (44) 

An 0(z)-correction is indicated, since £ Z/ g(d) from £u(d) produces such a term. In 
conclusion 

e = e (0) + z€ (1) = -(V <8> u + u ® V) + u>d 

- z Q(V <8> Vw + Vw ® V) - d 2 w - £ u (d)^ + 0(z 2 ) 
= -(V ® u + u <g> V) + wd 

- z (v ®Vw -d 2 w - £ u (d)) + 0(^ 2 ) . (45) 

This result agrees with that in jH] ( where d := — d ) when the Lie-derivative in (J45|) is 
written out in full. Using ^2] or performing similar calculations as in sect. 13.21 gives the 
standard result 

£ u (d) = (d ab . c u c + d c a u b , c + d c b u a , c ) 4> a <g> 4> h (46) 

with the covariant derivative corresponding to the mid-section. This Lie-derivative now 
measures the variation of the second fundamental tensor, the curvature tensor, with 
respect to membrane fields. It represents the change of curvature when only stretching. 

First versus second fundamental tensor To summarize, the Lie-derivative of the first 
fundamental tensor a ab produces parts of the leading in-plane strain e*- -*, whereas on 
the second fundamental tensor d ab parts of the bending tensor [HI E] introduce 

various measures of variation of the second fundamental tensor and linearize these. The 
same results are obtained here formally, where the Lie-derivative eventually acts on the 
curvature tensor in the expansion of the metric. 



How unique is the bending tensor? The discussion in this text considers the variation 
of d = d ab (fi a cg> 4> b under deformation. Since the Lie-derivative does not commute with 
raising and lowering indices, in fact £u(g) = 2e, the various measurements of bending 
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defined from the variation of d a b, d b a or d ab differ with terms of the form e 
instance [Oj considers the variation of d\. 



12 

d p. For 



What dominates: bending or stretching ? The most important terms in the strains ()45j) 
at short wave lengths, the principal part, can be identified as the highest derivatives. 
The term of the strain proportional to the vertical displacement has as its principal 
term V (g) Vw associated with bending and arises from the zf3 contribution to the 
in-plane strain. On the other hand, the principal part of comes from the tangential 
stretching field described by the V U, the in-plane strain. 

4. Shell dynamics 

In this section the stresses and moments in the shell are found. We follow [Oj but use 
exterior calculus instead. 

4-1. Force conservation in the bulk 

Forces acting on general surface elements dS are described in continuum mechanics 
|13| IT4*] in terms of the stress tensor <x. Thus the infinitesimal force on a surface element 
equals 



where X is an external volume body force. In this section we focus on the left-hand side 
of (|48[). that is the dynamic part of the equations of motion. 

In exterior calculus it is convenient to collect the forces on an infinitesimal element 
via the stress-form [TT] 



df = dS • a- . 



(47) 




(48) 



r = (ds ■ a-y 



a H dSi + c7 2i dS 2 + a 3i dS 3 . 



(49) 
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This form is a two-form and when integrated over a surface it gives the corresponding 
force (here in direction i, which by convention is put as the last index in the a m ). In 
continuum mechanics P is infinitesimal as in (|47|) but in exterior calculus the notation 
does not indicate so as this form is not exact in general. 
The total force is written as 

f = f i e i = f i <g>e i (50) 

Its covariant derivative is simple 

Vf = VP <g> e 4 + f ® Ve^ = VP ® e. t (51) 

as the frame and its dual are constant and similar for X, ii. 
In differential form ()48|) reads for the component % 

VP + XMV = pti'dV . (52) 
4-2. Covariant derivative of stress form 

First we discuss tangential directions. Hence, make a split in tangential and normal 
(indices a, b G {1, 2} versus 3) in (|52"j): 

VP = V(a la dS! + a 2a dS 2 + a 3a dS 3 ) (53) 
= d 1 (a la A 2 (l + zjR 2 ))da x A da 2 A dz + d 2 (a 2a A 1 (l + z/R^da 2 A dz A da 1 
+d 3 (a 3a (l + z/R 1 )(l + z/R 2 )A 1 A 2 )dz A da 1 A da 2 + w a 6 A P + u: a 3 A f 3 

Consider the first component a — 1. A simplification occurs in the last two terms of 
(|53[). since e.g. uj a % = ® a /(R a + z) restricts the relevant term in f 3 to <r a3 dS a by Q. 
Likewise the second last term has 6 = 2, where 

u\ = -1— (A ll2 eV(i + z/Rt) - A 2il e 2 /(i + 2 /i? 2 )) (54) 

A\A 2 

picks out corresponding terms of f 2 , that is a 21 for the first term in (|54|) and cr 22 for the 
second. 
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4-3. Momentum conservation 

The degree of freedom in the direction perpendicular to the shell is eliminated by 

rh/2 

)z=-h/T 



integrating ()52|) over the thickness coordinate, f^/ 2 



4.3.1. Integrated stress tensor Consider the force per length with respect to an edge 
along the middle section. For example, assume the edge is perpendicular to the first 
direction. That force density equals, 



f h/2 dSAz)^ rh/2 
N la = J ~ h/2 A - - = / dza la (l + z/R 2 ). (55) 

J-h/2 



-h/2 
'2 

A 2 da 2 

The notation is such that the first index denotes the direction to which the perpendicular 
edge is considered. Thus, N la is the force density along the edge 2. Likewise Q 2 = N 23 
is the force density perpendicular to the shell in the z-direction along edge 1. 
In the thickness integration produces precisely such terms: 



N a = N aa 

rh/2 

N la = / dza la (l + z/R 2 ) 



N 



2d 



h/2 
h/2 



/ dza 2a (1 + z/Rt) 

J-h/2 



Q 



-h/2 
-h/2 



1 = / dza xz (l + z/R 2 

J-h/2 



fh/2 

Q 2 = dza 23 (l +z/Ry) . (56) 

J-h/2 

referred to as stress resultants. Also corresponding moments are defined as: 

-h/2 



la 



M 



M 2a 



-h/2 

following the conventions of [3]. 

The resultants satisfy the constraint 



Mil A 

\ dza la z{\ + z/R 2 ) 

J-h/2 
rh/2 

/ dza 2a z{\ + z/R-y) (57) 

J-h/2 



N 12 -N 21 = -^M 21 --^M 12 . (51 
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(|58|) follows from the symmetry of the stress tensor 

= [ dz (a 12 - a 21 )(l + z/#i)(l + z/R 2 ) = N 12 + i- M 12 - N 21 - i- M 21 . (59) 
We express the resultants in tensor form as 

rh/2 

N ab = / dz (5 a c + z d a c ) a cb 



h/2 
h/2 



Q a = [ dz (5* + zd a c )a c3 

J -h/2 



-h/2 

rh/2 

M ab = / dz z (5 a c + z d a c ) a cb (60) 

J -h/2 

where a modified curvature tensor d is introduced as 

d = 4-0 1 ®0 1 + 4-^ 2 ® < / )2 = ( trd ) a - d - ( 61 ) 
tii til 

For other resultants, see [Hj. 

4-3.2. Force balance for shell element We continue the example from sect. 14.21 When 
integrating over z and multiplying with AiA 2 note that A a and R a are independent of 
z. Also the symmetry a 2 i = o~vi is used to get a term N 12 . The integration over the 
c?3(= d z ) -term is evaluated at the boundary. Collecting terms proportional to do 1 Ada 2 
gives the first equation of motion: 

d 1 {A 2 N 1 ) + d 2 (A 1 N 21 ) + A 1>2 N 12 - A 2A N 2 + A 1 A 2 (Q 1 /R 1 + q x ) = A x A 2 phu x (62) 
with the effective load 

rh/2 

= [ a 13 (l + z/RtXl + z/R 2 )] h J" + / dzX 1 (1 + z/R^l + z/R 2 ) . (63) 

1 J-h/2 

Similar calculations are done for the second component a = 2 in (J53|) . 
4-4- Moment conservation 

The moment in the direction z A is found by multipliying (|52jl with z. As before we 

consider i = a = 1 and the derivation is very similar to the above and differs only in a 

partial integration of the corresponding c^-term 
rh/2 

/ dzzd^a^il + z/R^il + z/R^A^da 1 Ada 2 (64) 

J-h/2 
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[za 31 (l + z/R x ){\ + z/R 2 )] h /*_ h/2 - j ' dza 31 (l + z/R 2 )j A 1 A 2 da 1 A da 2 



Ml 

- dzcx 31 (l + z/R 2 )z/R 1 AiA 2 da l Ada 2 

J-h/2 



-h/2 

The last term cancels with 



h/2 nh/2 

z uK A f 3 = / z uK A a 31 @ 2 Adz. 



3 A I = I 3 

z=-h/2 Jz=-h/2 

The second term in the parenthesis in (|64|) gives the integral for Q 1 , whereas the 
boundary term is grouped together with body moments as done above for the effective 
load q l . 

4-5. Normal resultants 

The final stress resultant involves the normal stress and is denoted Q a . The calculation 
can be done as in (|6*2j) and stated with derivatives of products of Lame-coefficients A, 
and resultants. 

However, the derivatives of the products can be expanded and written as covariant 
derivatives of the mid-section using ([28)1 . This can also be seen at an earlier stage in the 
calculation by keeping the mid-section frame {(f) 1 , </> 2 , dz] and evaluating derivatives 
using (125(1 . For simplicity we show this alternative approach for Q a but similar 
calculations can be done for the other resultants. 

The covariant derivative becomes 

Vf 3 = df 3 + u 3 a A P (65) 
with 

df 3 = d{a 13 @ 2 A dz) + d{a 23 dz A 1 ) + d^B 1 A @ 2 ) 

= (di(a 13 (l + z/R 2 )) + d 2 (a 23 (l + z/R^) + ^(^(l + z/R^l + z/R 2 )))ct> 1 A 4> 2 A dz 

-(a 13 (l + z/R 2 ))u\ A (f) 1 A dz + (or 23 (l + zjR x ))dz A u\ A 4? 

= [9i(a 13 (l + z/R 2 )) + ^(^(l + z/R x )) + 9 3 (a 33 (l + z/R x ){\ + z/R 2 )) 

+ (a 13 (l + z/R 2 )) f 2 21 + (a 23 (l + z/R x )) f\ 2 ] (fi 1 A <f) 2 A dz (66) 
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and 

u>\ Af = A a ba dS b = --l—a ba 5 ab dV (67) 

R a + z R a + z 

so 

J Vf 3 = (diQ 1 + f\ 2 Q 2 + d 2 Q 2 + r^g 1 - d ab N ab 

+ [a 33 (l + z/fliXl + z/i? 2 )]^*) 1 A <f> 2 = \v a Q a - al ab N ah + q 3 ] cp 1 A cp 2 . (68) 



z=-h/2 / 

as r*j = t^w] = for all i,j (no summation) by the antisymmetry of the frame 
connection coefficients. On the other hand, the kinematic body force equals 

rh/2 

/ pii 3 dV = phw (p 1 A <p 2 + 0{h 3 ) . (69) 

Jz=-h/2 

Thus to leading order 



V a Q a -d ab N ab + q 3 = phw. (70) 
4-6. Dynamical equations of motion 

In summary, using exterior calculus and moving frames one finds directly the shell 
equations in classical form [Sj: 

d 1 (A 2 N 1 ) + d 2 (A 1 N 21 ) + A 1>2 N 12 - A 2A N 2 + A 1 A 2 (Q 1 /R 1 + q 1 ) = A x A 2 phu x 
d x {A 2 N 12 ) + d 2 (A 1 N 2 ) + A 2jl N 21 - A li2 N 1 + A 1 A 2 (Q 2 /R 2 + q 2 ) = A x A 2 phu 2 
d 1 (A 2 Q 1 ) + d 2 (AxQ 2 ) - A 1 A 2 (N l /R 1 + N 2 /R 2 ) + A 1 A 2 q 3 = A x A 2 phib 

3 X {A 2 M X ) + d 2 (AiM 21 ) + A lt2 M 12 - A 2A M 2 - A 1 A 2 (Q 1 - m 1 ) =0 
d x {A 2 M 12 ) + d^M 2 ) + A 2A M 21 - A X)2 M l - A X A 2 {Q 2 - m 2 ) =0. (71) 

The right hand sides, the kinematic part, are also quantities integrated over thickness. 
They will be discussed in sect. 14.71 

Yet, the equations for the resultants can also be written covariantly: 

V b N ba + d a b Q b + q a =phu a 
V a Q a -N ab d ba + q 3 = phw 

V b M ba -Q a + m a = , (72) 
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since in ()71|) in the equations for N ab and M ab all the first four terms come from a 
tangential covariant derivative when we recognize the connection coefficients given in 
(|28|). Likewise for Q a as demonstrated in sect. 14.51 

Eliminating the normal stress Q a and assuming no body forces and moments gives 

V b N ba + d%V c M cb =phu a 

V a V b M ba -N ab d ba = phw. (73) 
4- 7. Inertia terms 

When the acceleration displacement field ii is integrated over thickness only those terms 
constant in z survive leading to the right hand sides of the three first equations in (fTTj) . 
But for the last two, the moment equations, (|52|) is already multiplied with z and only 
the acceleration of (3 l remains. This leads to terms 

which go as the cube of thickness and are neglected in this treatment. 
4-8. Plane stress 

Although all that we wanted to consider concerning moving frames and shells was 
presented in the previous sections we mention for completeness one final subject, which 
is how to relate the resultants with the displacements. Here a link between stress and 
strain is needed. These are in the form of constitutive equations and can be thought of 
as a generalization of Hooke's law. 

In the analysis of plates and shell the so-called plane stress approximation is used 
©• I n that approximation, the stress-tensor has almost the same form as in bulk 
elasticity except that the coefficients of elasticity are slightly altered. We refer the 
reader to IE] for a full discussion and just state the approximation in the case of 
isotropic materials 
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Following the notation of [TS] a 4-tensor is introduced 

H abcd = ( a ac a bd + a ad a bc^ + y fl ab flC d ^ 

such that 

ab zjabcd , t<7Cl\ 

v = z. z-H e cd . (76) 

1 — z/ z 

4-8.1. Stress resultants We are now prepared to integrate (|Hfi)l to get the resultants. 
Separate the stress tensor according to the different orders in the expansion in z: 

a ab ^ = -JL- H abcd (77) 

for % — 0, 1 with given by (|45)l . Plain calculation from (J56|) and (|71j) gives the 
resultants as 

N ab = ha ab(0) + ^ja cb(l) 

12 c 

M ab = ^ (a^ + d a c a cb ^) (78) 
using d defined by (|61|). Or, in terms of the plane stress tensor _f/" a6c,i 
jV ab = c # abcf efj + 5 H cbfe ef e 

M ab = B [H abcf e$ + d a c H cbfe ef e .) (79) 

The constants C and S are the stretching and bending rigidities 
Eh , „ £/i 3 

r- With B = ; 

1-z/ 2 12(1 -v 



C =^—2 with fi = 7^7— 3T- ( 8 °) 



In ()78j) the factors z/Ri from ()56j) are included leading to corrections proportional 
to h 3 /Ri, see jSlE!- If included the relation (|HH)l among the resultants is satisfied. Thus, 
some tensor algebra on (|78j) shows 

e ab N ab = -e ab d a c M cb . (81) 

Therefore the antisymmetric part of N equals that of d • M and in lines of curvature 
coordinates has the form 

N 12 -N 21 = -j- M 21 - i- M 12 (82) 
R2 R\ 

which coincides with (|58|) . 
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5. Summary and conclusion 
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The method of moving frames and exterior calculus allow a fast derivation of the 
equations of motion for a curved elastic shell. First, the kinematic stretching and 
bending field corresponds to a particular strain tensor field obtained using the Lie- 
derivative on the metric. Second, the equations describing the dynamics are found using 
stress forms which are differential two-forms encoding the stress tensor in a convenient 
way. Finally, all equations are seen to be covariant and hence valid in any coordinate 
system. 
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